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This article investigates an extension of General Relativity based upon a 
class of lifted metrics on the cotangent bundle of space-time. The dynam- 
ics of the theory is determined by a fixed section of the cotangent bundle, 
. representing the momentum of a fluid flow, and Einstein's equations for 

^ 1 the fluid applied to the induced space-time metric on the submanifold 

^ | of the cotangent bundle defined by the image of the section. This con- 

, struction is formally analogous to the extension of Galilean Relativity by 

f"^ ' Special Relativity, and is shown to reduce to General Relativity as the 

gravitational constant approaches zero. By examining the consequences 
of the model for homogeneous cosmologies, it is demonstrated that this 
construction globalizes the equivalence principle, in that, the perfect fluid 
^-l model of Special Relativity is sufficient to predict both the inflationary 

I and the current era. 
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1 Introduction 



In the standard formulation of General Relativity (GR), it is axiomatic that 
frames are non-dynamical objects determined by the choice of observers. This 
feature of GR is a consequence of the clock hypothesis axiom, which postulates 
that the rate at which an ideal clock (one that is corrected for the mechanical 
effects of acceleration) operates depends only on its 4-velocity. However, the 
role that frames play in dynamical structure of GR is more complicated than 
suggested by their place in the axiomatic framework. This is because matter, 
as it is modeled in the theory, often distinguishes a frame which becomes a 
dynamical object. Examples of such frames are the 4-velocity streamlines in 
fluid models and stationary observers in static space-times. The dependence of 
the space-time metric on the distinguished frame of the matter flow suggests an 
extension of GR in which the frame dependence of the metric is incorporated 
into the axiomatic structure of the theory. 

To see the possibility of such an extension, a frame of observers must be 
viewed as a geometric object in its own right. This is best accomplished if 
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a frame, rather than being viewed as a choice of coordinates on space-time, 
is represented in an abreviated form, as a closed 1-form field A on space-time 
M. The data that A carries concerning the associated field of observers is the 
nature of the observers' clocks, and the local decomposition of M into space- 
like submanifolds. Now a 1-form field can also be viewed as a section of the 
cotangent bundle T*M of M. From this perspective the geometric object most 
easily associated with a frame is its image A(M). 

To arrive at the proposed extension of GR, the frame of observers A(M) 
is thought of as an analogue of an observer's world-line in Galilean or special 
relativity. The fundamental problem in special or Galilean relativity is to fix a 
representation of time along a world-line. A similar question arises in extending 
GR which is to establish method for defining a space-time metric on A (M). 
The Galilean approach advocates an extrinsic procedure in which there is a 
fixed correspondence that takes A(M) to a model space where distances are 
measured. In Galilean relativity the fixed extrinsic structure is the universal 
time line R. In the present case, M and the bundle map tt: T*M —* M serves 
this purpose. On the other hand, relativistic approach postulates that distance 
measurement is intrinsic to the submanifold A (AT) itself, and may depend on 
the position of the submanifold in the ambient space. 

The Galilean approach can be seen to yield the standard structure of GR. It 
is encoded in the axiom which for emphasis Einstein states as continuum spatii 
et temporis est absolutum |Q. In this article we shall investigate whether the 
relativistic approach may be the more powerful. Following Einstein's dictum, if 
the relativistic approach is to be a viable alternative, it must be demonstrated 
to be consistent, in the appropriate limit reduce to prior tested theories, and 
have interesting consequences. Consistency is of course a very difficult question 
in physics. The second criteria of reduction to existing theoretical structure at 
least suggests a level of consistency similar to existing theories. In following we 
hope to demonstrate that an extension of GR using the relativistic approach to 
measuring space-time distances satisfies the last two of these criteria. 

To show that the relativistic approach does in fact have interesting conse- 
quences, we shall examine the implications of the extended theory for homoge- 
neous cosmological models. Here, because in the extended theory the metric is 
frame dependent, it is most natural to apply the equivalence principle in the 
matter frame. This leads to a third rather than second order equation for the 
cosmic scale factor; the solutions of which appear to have interesting global con- 
sequences. It will be shown that if the universe is assumed to evolve towards a 
dust state, then it must have originated from an inflationary state. In fact the 
model suggests that the alternation of standard and inflationary phases may 
be a cyclic process. Analysis of the model shows that it has the same conse- 
quences as current homogeneous inflation models but with the advantage that 
cosmic evolution from the inflationary phase to the present can be described 
as a single process. The mathematical structure presented here does not yet 
apply to the study of fluctuation in initial cosmic matter density, and so it can 
not be expected to account adequately for phenomena such as reheating and 
condensation. 
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2 Extended General Relativity 

In this discussion a frame on a space-time M will be represented by the pair 
(A, V), where A is 1-form field on M, and V is vector field on M. A frame 
provides a decomposition of the tangent space to M at p that is given by 
TM p = span(V p ) (J) ker(A p ). This decomposition is interpreted as the infinites- 
imal decomposition of TM p into space and time directions by the observer as- 
sociated with the frame and located at p. Note that this decomposition is 
independent of the choice of scale for the vectors V and A. The scale of V and A 
has useful descriptive properties. The scale V is associated with the observer's 
choice of time coordinates and the scale of A can be interpreted as a mass scale 
for particles. 

Although the association of 1-forms with mass is a standard correspondence 
in GR that is usually established from momentum consideration, we provide here 
an argument for this relation based on Newton's Second Law. Suppose that W 
is another vector field on M, then the invariant velocity w of W in the frame 
(A, V) is given by w = (1/X(W))W. Note that w is invariant quantity in that 
it is independent of the choice of scale of W, and hence is a geometric objects 
associated with the flow lines of W. The velocity w, however, is dependent on 
the choice of scale of A. Newton's Law in the given frame may be expressed as 
mw — F(w), where F is the force and w is the invariant derivative of w. The 
invariant derivative is given in terms of the Levi-Civita connection V by the 
expression 



Here ± is the projection onto the complement of V. 

Now suppose that the frame (A, V) is replaced by (kX, V) where k is a positive 
constant. Let w' be the invariant velocity of W in the frame (k\,V). The 
above definitions imply that mw' = m{l/k 2 )w = (l/k 2 )F(w) = (l/k)F(w') = 
or kmw' — F(w'). Hence in the frame (k\,V) the mass of W is m! = km. 
Hence A can be imputed to have the same dimension as m 

The equivalence principle guarantees that relativistic dynamics can be re- 
duced to Newtonian dynamics infinitesimally along any flow line of V. However, 
it is only under very strong conditions that this correspondence can be extended 
locally. In the case where A = m£(V), £(V) being the 1-form metrically dual to 
V, the Newtonian constructions gives an equivalent description of metric geom- 
etry only when V is a gradient Killing field Q . This assumption implies that V 
is holonomically trivial or R(-, -)V — 0. Unfortunately, a self gravitating frame 
of massive observers cannot be holonomically trivial, and this fact leads to a 
global incompleteness in the description of the physics of such systems. Math- 
ematically, this difficulty surfaces in the lack of a globally holonomic frame in 
which to apply the equivalence principle. This mathematical fact is mirrored 
in physical principle that gravity destroys states of local thermodynamic equi- 
librium by invoking the other fundamental forces and thus altering the local 
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equation of state. This is particularly true in cosmology where the evolution of 
fluid streamlines can only be determined by fitting equations of state appropri- 
ate to each of the local temperature regime found in the evolving universe. The 
fact that there is no universal law that describes cosmic evolution thus can be 
traced to the fact that we lack a universal description of the forces in nature. 
However, the fact that such a simple model as the perfect fluid model appears 
to describe the bulk of cosmological space-time, suggests there should exist an 
underlying principle from which a globally significant cosmological model can 
be derived. 

Here we shall investigate the possibility that such a principle may be derived 
from a generalization of the equivalence principle that arises in an extension of 
GR in which the space-time metric is frame dependent. This extension of GR 
is obtained by a procedure that is analogous to the procedure used to obtain 
special relativity (SR) from Galilean relativity (GLR). Recall that the Galilean 
theory of time can be formulated by supposing that the space-time manifold 
M possesses a fibration t: M — > R. This fibration is assumed to possess a 
connection represented by a vector field V on M which must necessarily have 
the property that if <f) s is the flow of V then (p Sl (t~ 1 (s 2 )) — t~ 1 (s\ + s 2 ). Thus 
the pair {dt, V) determines a synchronous field of observers on M with the 
special property dt(V) = 1. Let N be a model for the fiber of t. GLR sup- 
poses that N and R are separately endowed with the complete Riemannian 
metric tensors q and g, and that g is solely responsible for determining the 
time between two events p,q € M. If 5 g is the topological metric on R de- 
termined by g, then according to GLR, the time between p and q is given by 
6g(t(p),t(q)). In SR time does not have this well defined character but depends 
on the observer with respect to which these events are referenced. The metric 
that an observer uses is induced on the observer world line from a metric g on 
M that is obtained by summing g and q. This is done using the decomposition 
TM p — span(V^) 0ker(dt p ), the natural isomorphisms t* : span(V) — > TR, and 
the identification i s : t~ 1 (s) — > N. The metric g is then given by the expression 
g = —t*g@(l/c) 2 i*q where c is a universal constant. In the limit as c — ► oo, 
SR reduces to GLR. 

This procedure can be carried out again with the base manifold R replaced 
by the space-time manifold M, and the fibration t: M — > R replaced by a 
fibration ir: M — » M. The manifold M that replaces M can be discovered by 
noting that in SR the fiber t~ 1 (s) parameterizes the positions of all possible 
observers at time s. Thus, by analogy it is natural to suppose that for an event 
p e M the fiber n~ 1 (p) parameterizes the frames of all possible observers at p. 
The above description of frames indicates that the fiber should be TM p x T*M p . 
However, it will be assumed that there is a relation I: TM — ► T*M, and that 
the frames under consideration all have the form (A,^ _1 (A)), and thus the fiber 
of 7r at p may be assumed to be T*M p , making M = T*M. 

As in SR, a metric can be constructed on T*M using the Lorentzian metric 
on M and the induced metric on the fibers of T*M. To interpret this construc- 
tion it is helpful to recall how relativity unifies space and time dimensions. In 
the case of SR, tangent vectors to N are assigned the dimension of length L 
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while tangent vectors to R are assigned the dimension of time T. Metric ten- 
sors are supposed to be dimensionless. If the universal constant c is supposed 
to have the dimension L/T, then for vector fields V and W on M, g(V, W) has 
dimension T 2 . This suggests that the dimensional structure of M should be 
simplified by measuring the length of tangent vectors to N in T-units with 1 
L-unit equaling 1/c T-units. In such a system c = 1 and all physical tangent 
vectors to M have dimension T. 

In the case of T*M, the above discussion of Newton's Law suggests that the 
fibers of T*M should be given the dimension mass M. Thus in this case we see 
that there is the need of a physical constant 1/G of dimension M/T relating 
the base and fiber dimensions. The constant G is easily seen to have the same 
dimension as the gravitational constant. 

To construct a metric on T*M denote the space-time metric by g and pro- 
ceed analogously to construct a metric g on T*M. In analogy with GLR, the 
horizontal distribution H of the Levi-Civita connection plays the role of the 
field V. It is a well known fact in connection theory that TT*M decomposes 
as TT*M p = H p ®VT*M p , where VT*M is the vertical distribution of the 
fibration it: T*M — > M. Because T*M is a vector bundle there is a natural 
isomorphism i p : V p — > T*M n ( p ) defined for each p € T*M. The isomorphism 
i p together with the isomorphism 7r* p : H s — ► TM„r p \ can be used to define a 
class of metrics on TT*M with the property that H and VT*M are orthogonal 
distributions. Of particular interest is the metric g given by 

This expression is considerably more complicated than the corresponding 
expression in SR due to the presence of the metric dependent conformal factors. 
There are several reasons why this particular form was chosen for g. First, the 
parallelism of this metric can be shown to be closely related to Fermi parallelism 
along curves on the space-time Thus, the relativistic effects of acceleration 
are encoded in the geometric structure defined by g || . Second, since g(X, Y) 
should have a well defined dimension, the dimensions of the vertical and hor- 
izontal summand of this quantity should be equal. This implies that e/g(p,p) 
has dimension of the gravitational constant T/M, and so e has the dimension 
of Plank's constant MT. Thus the dimension of the conformal factors can be 
factored as MT/M 2 . This factorization is similar to the factorization of the 
gravitational constant G = h/Mp 2 where Mp is Planck's mass. This particular 
way of expressing the dimensional dependence of g will prove very useful in 
cosmological applications. 

To interpret the metric <?, consider a field of observers associated with a 
section A of T*M . Just as in SR where the world-line of an observer inherits 
a metric that measures time displacement, in this case the image of space-time 
X(M) inherits a metric G'glxiM) which measures the space-time displacement 
in that frame. Here a dimensional constant G is introduced with dimension 
T/M so that when evaluated Gg has dimension T 2 . Thus in extended general 
relativity (EGR), an observer field A detects not the space-time metric ?r*g but 
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rather the space-time metric Gfg\\(M)- It will be demonstrated in section 3 that 
if G — ► and e — > in such a way that Gm 2 /e — > 1, then EGR reduces to GR. 

3 Global Equivalence in Friedmann Geometry 

This section shall examine some consequences of EGR in the standard cosmolog- 
ical model. In this case, space-time M is modeled by the product RxJV where 
N is now a simply connected 3-dimensional space-form of constant sectional 
curvature k. If q is the constant curvature metric on N and if t is the natural 
affine coordinate on R, then consider the class of Friedmann metric tensors 

g = -dt 2 @S(tfq, 

where S(t) is assumed to be a smooth positive function. Frcidmann geometries 
possess a distinguished frame determined by the unit vector field T in the distri- 
bution TR x 0. This frame is described by the pair (T, s) where s = m£(T) and 
m is some intrinsic mass. To study the geometry of the frame (T, s) first observe 
that g\ s (M) is °f Robinson- Walker type. To determine the form of g on s(M) 
use the induced map s* : TM — ► Ts(M) to parameterize Ts(M). It follows from 
the theory of linear connections that for any tangent vector u E TM p , 

where V is the dual Levi-Civita connection of g and denotes the horizontal lift 
to TT*M determined by V. From the identities 

V T s = 

V x s = mlogS(t)£(X), 

where X is a vector field in the distribution x TN, and the above expressions 
for s»m and g, it is easy to see that 

g(s*T, s*T) 
g(s*X, s*Y) 

These expressions imply that since tt: s(M) — > R x N is a diffcomorphism 

g\s(M) = ^(n*dt 2 -(1 - 4,(log S(t)f)S(t)\* q ). 

Since the intrinsic geometry of the frame determined by s is a Friedmann 
geometry it provides a model of a self gravitating perfect fluid. To complete 
the model one needs in addition to Einstein's equations an equation of state 
for the fluid. At this point it is customary in GR to appeal to the equivalence 
principle and to SR for an equation of state. In EGR, since the metric is frame 



m 2 
e 

Ttl 2 f 2 

= (1-— {\ogS(t)f)S(t) 2 q{X,Y). 
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dependent, there is the question of choosing the frame in which to apply the 
equivalence principle. The hypothesis that shall be advanced is that the equiv- 
alence principle should be applied in the frame of the matter flow, and that in 
this frame it yields an equation of state with global validity. In other words, the 
entropy generating effects of gravity can be eliminated by "viewing" space-time 
from the matter frame. This idea can be thought of as an extension of Einstein's 
original version of the equivalence principle in which gravitational forces were 
eliminated by accelerated frames. The original form of the equivalence princi- 
ple was too strong as it only applies infinitesimally or to pseudo-gravitational 
forces as note above. However, because the metric in EGR is frame dependent 
the original form of the equivalence is no longer limited to flat geometries and 
appears to have interesting implications. 

To derive the consequences of this assumption recall that for a simple gas, 
kinematic analysis implies that the energy e, pressure p and density p satisfy 
the equations of state e — p = ap and that the energy momentum tensor has 
the form T = (e + p)V ® V + pg. When coupled with Einstein's equations, 
the equation of state yields a system of equations that determine S(t). For 
simplicity let S'{t) 2 = (1 - (e 2 /m 2 )(log'S(t)) 2 )S(t) 2 . From Einstein's equations 
in the frame determined by s, we have 

e = 1(008^))- + ^) 

V = -^(21ogS'(t) + 3(logS'(t)) 2 + ^). 

Further the matter conservation law implies that there is a constant g such that 
P = g/S'{tf. 

To obtain the equations for S(t) implied by these relations and the equation 
of state, first assume that, as suggested in the last section, G = e/m 2 . It will 
also be useful to transform to the universal time scale by setting u(t) = S(kt) 
where k = e/m. With these assumptions a calculation shows that u(t) satisfies 

2 

logw(l — logw )logu = 

i " n i ' 2 \i-\ 3 + 4a • 2 •• 2 3 + a ■ 2 

logutl — logu )(1 logw j — log it (1 — logu j + 

a 2a 

3(1 + 0), -2 • 2 o 3 + a nk 2 -2 1 k 2 Gq , • 2 1 

^-^log M (1-Iogu ) +—^T( 1 - l0 S u ) + ^ e — ( 1 - l0 S u ) 2 - 

It is physically reasonable to assume that the constants (3 + a/2a)nk 2 and 
(l/2a)k 2 Gg are both extremely small, and so the last two terms in the above 
expression can be neglected. If these terms are dropped then the above equation 
for u can be rewritten as a planar system in the variables x(t) — logu(i) and 
y(t) = logw(t) which after rescaling time has the form x = F x (x,y) and y = 
F y (x, y) where 

F x = x{\ - x 2 )y 

„ 2 \/i 3 + 4a 2 \ 2/1 3 + a 2 3(1 + a) 9/1 2 \2 
F y - y(l-x 2 )(l — x 2 )-y 2 (l-—— x 2 )+ v ^ > x 2 {\-x 2 ) 2 . 
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This system, although it can not be explicitly solved, can be analyzed using 
the qualitative theory of planar systems. To interpret this analysis it should be 
kept in mind that the xy-phase plane is related through Einstein's equations to 
the Hubble constant-pressure plane of the cosmological model. First note that 
the above system possesses four stationary points (0, 0), (—1, 0), (1, 0) and (0, 1). 
The stationary point (0, 1) is easily seen to be of hyperbolic type. However the 
analysis of this equation is complicated by the fact that (0, 0) is semi-degenerate 
and (—1,0) and (0, —1) are totally degenerate. The stationary point (0,0) will 
be shown to be topologically of hyperbolic type and the stationary point at (0, 1) 
will be shown to possess elliptic and parabolic sectors. Using the symmetry of 
the equation it may be assumed that x > 0. This is equivalent to assuming that 
the Hubble constant is positive. 

To understand the global behavior of solutions to this system first consider 
the structure of solutions near (0,0). Since F x (0,y) = 0, it follows that x — 
defines a solution curve that adheres to (0,0). To find the other adherent 
solution invoke polar coordinates by the standard relations r = F r {r, 9) and r6 = 
F e (r,6) where F r (x,y) = (xF x (x,y) +yF y (x,y))/r a,nd F g (x,y) = (xF y (x,y) - 
yF x (x, y))/r. The other adherent solution can be located by considering a sector 
bounded by y = and a solution of F r (x,y) = 0. Note that the equation 
F r (x,y) — is quadratic in y and may be solved to give radial isoclines 

(1 - X 2 ){1 - 2±^X 2 ) ± ^/(l - a;2) a; 4(3H^a _ 12 1±£L X 2) + (1 - x 2 ) 2 (l - X 2 ) 2 
V±= 2(1 -^X 2 ) 

We may distinguish radial isoclines by the property that y + (0) = 1 and y_ (0) = 
0. Let I be the graph of y- for x > 0. Consider a small sector W bounded 
by y — and Iq. The vector field F exits both radial boundaries of W. First, 
because F x (x, 0) = and F y (x, 0) > 0, F exits y = 0. Second, a Taylor 
expansion for Fg(x,y) developed along Io implies that Fg(x, y-(x)) = —x 3 + 
o(x 4 ), and consequently, if the sector W is chosen sufficiently small, Fg < on 
the segment of Io bounding W. As a result, there must either exist a unique 
solution or an interval of solutions that tend to (0, 0) as t — > oo with initial 
conditions on the arc bounding W. 

To see that there is a unique solution, examine the first order ode r(d9/dr) = 
^(x,y) = Fg(x,y)/F r (x,y). Introduce the angular operator 9 = (xd/dy — 
yd/dx)/r. It follows from a generalization of a lemma of Lonn that there is 
a unique solution adhering to (0,0), if 0\& < in the interior of W ||. To 
verify this inequality, note that 8^(a;,y) = Q(x, y) / F r (x, y) 2 , where Q{x.y) is 
a polynomial of order 14 in x and y with properties that Q(x, y) = —4ay 2 (x 2 + 
y 2 ) + o(r 4 ) and Q(x, 0) = -3(5a + 3)(a + 3)x 6 + o(x 7 ). These facts imply that 
Q(x,y) is negative definite in a neighborhood of the origin, which implies that 
if W is chosen small enough, then 9^ < on the interior of W . Hence, there is 
a unique solution that adheres to (0, 0) in W , and determines a stable manifold 
Jo for the stationary point (0, 0). 

To understand the global backward time evolution of Jo, first consider the 
isoclines F y {x,y) — 0. This equation is again quadratic in y and has the solu- 
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tions 

(1 _ 3±^x 2 ) ± ^\2^x 4 + {l-x 2 ) 2 \ 
2(1 -^ 2 ) J ■ 

These isoclines again can be distinguished by y~(0) = and y+(0) = 1. The 
most crucial feature of the graphs of y_ and y + is the position and sign of 
the vertical asymptotes. It can be seen that only y_ possesses an asymptote 
at x — y/2ct/(3 + a) while y + is finite at x and has the value y+(x a ) = 
(3(3 — a)(l + a))/(4(3 + a)(7a + 3)). The sign of the asymptote depends on 
whether xq > 1 or a;o < 1. It is easily seen that lim - y~(x) = — oo if xq < 1 
and lim x ^ x - j/-(x) = oo if Xo > 1. Also note that both isoclines have the 
property that y±(±l) = 0. 



y± = (1-x 2 ) 




Figure 1. Flow lines for the equation x = yx(l — x 2 ), y — y(l — x 2 )(l — Ax 2 ) — 
y 2 (l — (l/2)x 2 ) + (3/2)x 2 (l — x 2 ) 2 showing stable manifolds Jo and Ji and radial 
and cartesian isoclines Iq and 7i. 



To understand the relation between the stable manifold Jo and the isocline 
Ii determined by the graph of y_, we develop a Taylor expansion for Jo at (0, 0). 
We may assume that Jo can be parameterized as the graph of a smooth function 
f{x); where / must satisfy F y (x,f(x)) = f'(x)F x (x, f(x)). To obtain the first 
two terms in the Taylor expansion, assume that f(x) = Cx 2 + Dx 4 + o(x 5 ) 
and substitute into the above expression. It can be seen that f(x) = — (3(1 + 
a))/(2a)x 2 + (3(1 + a)) 2 / (2a) 2 x 4 + o(x 5 ) while a Taylor expansion for y_ at 
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(0,0) has the form y_(x) = -(3(1 + a))/(2a)x 2 - (3(1 + a)) 2 /(2a) 2 x 4 + o(x 5 ). 
Thus the stable manifold Jo lies above isocline I\ . 

Consequently, the global behavior of J is determined by the position of the 
asymptote xo- If a < 3 or xo < 1, Jo lies in an unbounded region defined by 
y = 0, x = 1 and Iq in which case F y > and F x > 0. Since x = 1 is a solution 
curve Jo cannot intersect it , and thus is asymptotic to it as t — > — oo. The case 
where a > 3 or xq > 1 is quite different. Here Jo lies in a bounded region U 
defined by y = and ii (see Fig. 1). In this case divide I\ into two segments 
I[ and I" at the point (xi,yi) where I\ attains its minimum y- value. Assume 
(0, 0) G I[ . It is easy to see that the vector field F exits U along y = and I[ 
and enters along I". Consequently, in backward time Jo either adheres to the 
stationary point (1,0) or leaves U along I". In latter case since F y (x,yi) < 
for xi < x < 1, and since x = 1 is a solution, it follows that, if Jo enters the 
region bounded by I", x = 1 and y = yi, it must adhere to (1, 0). Thus in the 
case where a > 3 the manifold J is an unstable manifold for (1, 0) and a stable 
manifold for (0,0). 

The unstable manifold J\ leaving the stationary point (0, 1) can also be seen 
to adhere to (1,0). This follows from the fact that the boundary of the unit 
square [0, 1] x [0, 1] cither consists of solutions curves, x = and x = 1, or 
segments where the vector field is entering. Since J\ enters the unit square, it 
must adhere to a singular point on the boundary. However, it cannot adhere to 
(0, 0), since the only inward pointing characteristic direction with angle between 
and 7r/2 is uniquely associated with J . Thus Ji must adhere to (1, 0). 

The final step in understanding the global behavior of the flow of F is to 
analyze the stationary point (1,0). Here we shall restrict our attention to the 
case where a — > oo; that is to zero pressure cosmologies. In this limit we find 
that 

F x = x(l - x 2 )y 

F y = ? y(l-x 2 )(l-4x 2 )- y 2 (l-ix 2 ) + ^x 2 (l-x 2 ) 2 . 

To study the stationary point at (1,0) introduce the coordinate change 
z = 1 — x and transform to polar coordinates letting z = rcos(9) and y = 
rsin(#). After the standard rescaling of time, one obtains the equation for the 
blow-up of the vector F at (1,0); namely f = R(r,9) and 9 = S(r,9) where 
R{r, 9) = rR (9) + n(r, 9) and S(r, 9) = S (9) + £(r, 9) with lim r ^ £(r, 9) = 
and lim r ^o v(' r i 0)/ r = 0- The characteristic directions are determined by the 
equation So(0) = 0. The solutions to this equation are easily obtained since 

15 

S (9) = — cos(0)(l +cos(26> + </?)) 

1 K 7 

Ro(9) = —sm(9)(—+cos(29 + <p)) 

where ip = arctan(4/3). The equation So(6) — is easily seen to have solutions 
7r/2, — 7r/2, (ir—Lp)/2 and (ir— f)/2+TT. Also a short calculation gives (ir—<p)/2 = 
arctan(2). 
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To determine whether there are solutions to the ode adhering to (1,0) along 
these characteristic directions, first note that at any of these angles R(&) ^ 0. 
Consequently, the existence of such solutions is determined by the Taylor expan- 
sion of So at the given characteristic angle. In general, if the first non- vanishing 
coefficient is of odd degree then there exist solutions to the ode adhering along 
that characteristic direction. However, if the first non-vanishing coefficient is 
even then there exist either no solutions or an interval of solutions adhering 
along that characteristic direction. Examination of the above expression for So 
shows that the Taylor expansions of So at n/ 2 and — ir/2 are of odd degree, and 
at arctan(2) and arctan(2) + it they are of even degree. In both cases there are 
solutions adhering to (1,0) along these characteristic directions 0. 

Any solution with initial conditions (x,0) for < x < 1 must adhere in 
forward time to (1, 0) along the characteristic direction arctan(2) or 7r/2, and in 
backward time along the characteristic direction —tt/2. Since solutions adhere 
to arctan(2) exist, there must exist a closed interval [xo, 1] such that if x € [xo, 1] 
then solutions with initial conditions (x, 0) adhere to (1, 0) along the arctan(2) 
direction. Numerical results suggest that in fact x\ — 0, and that the last 
solution adhering to (1,0) at the angle arctan(2) is J\. Hence it is conjectured 
that the solution curves x = for < y < 1, Jo and J\ form a separatrix for an 
elliptic sector of the stationary point (1, 0). 

To develop a physical interpretation for these results, first note that the 
rescaling u(t) = S(kt) and Einstein's equations imply that points in xy-plane 
are related to the energy and pressure of cosmological states by the relations 



The most important integral curve of the vector field F for the interpretation 
of these results is J . It not only approximates the standard cosmological model 
for large time, but also, in the case where a > 3, it originates from an inflationary 
state. To see this, rescale vector field F so that its parameter correctly reflects 
the time evolution of the original problem; that is x = y. This rescaling can 
be accomplished by representing the solution curve as the graph of f(x), and 
then solving for x(t) by integrating x = f(x). In the case of Jo at (0,0), from 
above, we know that / has the form f(x) = — ((3 + 3a)/2a)x 2 + g(x) where 
\im x ^o g(x)/x 3 = 0. Consequently, integration gives that x(t) = (2a/ (3 + 
3a))(l/t) + b(t) where limt_>oo b(t) = 0. But, this is just the asymptotic behavior 
of solutions to the standard model with sectional curvature k = and density 
Thus Jo or solutions close to Jo are candidates for equations of state that 
give reasonable cosmologies. Following J in backward time, it either originates 
from a state of infinite pressure in the case where a < 3, or a state where 
p = — e in the case where a > 3. The latter equation of state is the equation for 
inflationary cosmology. The fact that for a > 3, cosmological models describe 
cold matter weakly interacting with radiation indicates that, in this scenario, 
the inflationary phase is predicted by the current condition of the universe. Note 



e 




P 



-^(2y + 3x 2 ). 
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that the energy density at the singularity is 3m 4 /e 2 . If m is taken to be the 
Planck mass this is consistent with the chaotic inflation hypothesis. 

To gain a better insight into the significance of this model for inflationary 
cosmology, requires a more detailed examination of the solutions near the sta- 
tionary point (1,0). Using a similar analysis as was applied at (0,0), it can 
be shown that solution curves that adhere to (1,0) are the graphs of either 
f(z) = —z 1 / i C{z) where C is analytic and C(0) > or g(z) = 2z + h(z) where 
lim z _^ h(z)/z = 0. These solutions correspond to the — n/2 and arctan(2) char- 
acteristic directions respectively. To find x(t) solve z = —f(z) and set x = 1 — z. 
In first case with the initial condition (to, zq), integration gives for < z < zq, 
z 3 / 4 C(z) - z 3 /4 C{z ) = t - 1 whcrc^C(O) = 4/(3C(0)). Let u 4 / 3 D(u 4 / 3 ) be the 
inverse of z 3 / 4 C{z) with D(0) = 1/C(0) 4 / 3 , then the solution x(t) has the form 

x(t) = 1 - ((t - t ) + zJc(zo))*I>(((t - t ) + zlc(z )) 4/3 ). 

One significant observation that follows from this expression is that, if Jo is pa- 
rameterized by the natural time parameter, then the stationary point is reached 

3 ~ 

along Jo in finite time. In fact, since z 4 C(zq)) > the singularity is reached at 
time t\ < t , t\ = to — z^ 4 C(z ). Knowing x(t) we may also calculate the scale 
factor u(t) = x(t)/x(t). Integration gives 

u(t) = u(tx) exp((t - h) - ((t - t ) + zlc(z ))iD(((t - t ) + z|c(z )) 4 / 3 ). 

The analysis of this expression shows that, although for small t — ti, u(t) grows 
exponentially, exponential growth does not persist long enough to account for 
inflation. To find true inflationary growth in this model we examine the behavior 
of solutions adhering to (1,0) along the characteristic direction arctan(2). 

For solutions adhering along the arctan(2) direction, the integral of z = 
—g{z) can be expressed in terms of the initial data (to,zo) and a bounded 
function b(z), and has the form zb(z) = zob(zo)e~ 2 ( t ~ to * > for < z < z . If the 
inverse of zb(z) is vd(v), where d is also bounded, then x(t) has the form 

x(t) = 1 - z b(z )e-^ t - ta U(z a b(z )e- 2 ^ t «y). 

Hence the scale factor u(t) has the form 

u(t) = exp((t - t ) + z a b(z )e- 2 ^d(z b(z )e-^ t - t °y)). 

From this expression we see that as t — > oo the scale factor undergoes exponen- 
tial growth. This fact suggests that the solutions of the elliptic sector may be 
given physical interpretation if it is supposed that the inflationary phase of the 
universe's development was the final state of a previous cycle. 

As final remarks, observe that both the flatness of space and the adherence 
to the critical density can be deduced from this model, but in a manner quite 
different from the way in which these conclusions arise in standard inflationary 
models. First, the fact that the density of matter must be extremely close to 
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critical density is a consequence of the age of the universe relative to universal 
time scale k. Since, if cosmic evolution did not follow a path extremely close 
to Jo, it would have suffered a catastrophe before its current age. The flatness 
of space is a consequence of the same scaling law, because for large time scales 
the curvature may be neglected as it is multiplied by a factor of k 2 . Finally, 
observe that the stationary point at (1, 0) corresponds to singularity in g\ s tM)- 
Notice that if x 2 = 1, then S'(t) = 0. At such points Ts{M) is tangent to the 
null cone of g and so space-time at the singularity does not have a well defined 
spatial metric. 



4 Submanifold geometry for a geodesic observer 
field 

In the previous section the implications of EGR arc studied for the simplest 
interesting example, namely Fricdmann geometries. In general it is difficult 
to compute the intrinsic geometry of an arbitrary submanifold s(M) of T*M. 
Although general formulas for the curvature of such submanifolds can be de- 
veloped, they are not at this point particularly illuminating and so will not be 
fully exposed here. Rather, this section will demonstrate that if s has constant 
length and i~ 1 (s) is a geodesic vector field, then the difference between the 
curvature on s(N) induced by ir*g and g depends only on the first three co- 
variant derivatives of s, and can be developed in a power series in e with first 
non- vanishing term of order e 2 . Thus, the additional terms that appear in the 
Einstein tensor when g is used in the place of it* g can be viewed as additions to 
energy momentum tensor that are coupled to gravity by higher powers of the 
gravitational constant. 

To understand the intrinsic curvature of s(M), first calculate the Levi- 
Civita connection V of g on T*M. Let R be the curvature tensor of the 
dual Levi-Civita connection in T*M . By metric duality R determines a ten- 
sor R defined on T*M x T*M x TM with values in TM by the expression 
g(R(s,t)v,w) = g(R(v,w)s,t). For vector fields X and Y and 1-form fields s 
and t, V is determined at a point p 6 T*M by the expressions 

V^F, = v7Y+-g(p,p)g(X,Y)i- 1 P -i- 1 ^R(X,Y)p 



g(p,p) 2g(p,p) 2 



Vi-i t i _1 s p = - — {g(p, f + g(p, s)i~H - g(s, tji^p). 



Now fix a section s of T*M that has constant length and is closed. Since V.s 
is a bundle morphism between TM and T*M, define V.s T to be its transpose; 
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that is, V.s T is a morphism between T*M and TM. The conditions that s be 
of constant length and closed can be stated in the form V s .s T = and V[/s = 
(S7ius T where U is a vector field on M and £: TM — > T*M is the identification 
map determined by the metric g. Note that these identities immediately imply 
the classical result that i~ x s is a geodesic vector field. If s is assumed to be 
time-like and to 2 = —g{s, s) then for vector fields U and W on M the above 
expressions for V imply that 

v s ,c/s*w = v d+i _ 1Vus w + r 1 v w s 

= VuW + ^ (R(s^u)W + R{s^ w )u) - ^i-'RiU, W)s + 
/to 2 1 \ 

rVt/VffS - ( -^g{u, w) + -^g{^ua, v w s)\ r 1 *. 

After some calculation using this form for V s ^uS*W, one obtains the following 
expression for the sectional curvature K(s*U, s^W) — g(R(s*U, s*W)s*U, s„W) 
of g along s(N) 

K(s*U,s*W) = 

9 

TYl 6 / 

K(U, W) H j ( 9{VuR(U, W)s, V w s) + 

e to V 

<?(V w i?(W, C/) s , Vc/s) + s(i2(I7, W)V v a, V w s)) + 

-^9(R(U, W)a, R(U, W)s) + (g(R(s, V uS )W, R(s, V w s)U) + 
g(R{s, V u s)W 1 R(s, W v s)W) + g{R(s, V w s)U, R{s, V w s)U) - 
3g{R{s,V u s)U,R{s,V w s)W)) + (g{V uS , V w sf - 

/ TO* V 

g(Vua, V uS )g(Vws, V w «)) - ^ (g(U, W) 2 - g(U, U)g(W, W)) . 

To obtain the intrinsic sectional curvature of s(M) from this expression 
requires the second fundamental form of s(M), a(U,W) = P /± (V s *uS*W). The 
orthogonal projection P f± on to the compliment of Ts(M) can be expressed in 
terms of covariant derivatives of s. To see this note that the tangent space to 
s(N) and its compliment are described as the graphs of the endomorphisms 
C: H s ^ VT*M S and C* : VT*M S -> H s ; that is, Ts(M) s = H S ® C{H S ) and 
Ts(M)j = VT*M S C(VT*M S ). It is easy_to see that for a vector field U, 
C{U) = i^VuS and C*{i-H) = (e 2 /TO 4 )V t s T . The projection P'^ can be 
expressed in terms of C and C* , and the orthogonal projections P and P" 1 that 
map TT*M onto H and its orthogonal compliment VT*M, and has the form 

p/-L = ( P ± + c*p^)(p^ _ CC^P- 1 )*^- 1 - CP) 
where denotes the group inverse. 



G. Martin Inflation and Global Equivalence 



15 



The assumption that s is of constant length and closed implies that £V.s has 
the structure of a 3-dimensional symmetric linear map A, and hence satisfies 
the polynomial identity A 3 = ^\A 2 — ^ 2 A + 6- The fundamental invariants 6, 
£2 and 6 are given in terms of the traces of powers of A by the expressions 

6 = tr(A) 

6 = i(tr(.4) 2 -tr(A 2 )) 
6 - itr(A) 3 - ltr(yl)tr(A 2 ) + itr(A) 3 . 
These invariants may be used to reduce (1 — A 2 )^ 1 to the expression 



l+e 2 


-26 


+66 + ^16- 

6 - 66 


|666 + 


3^3 




-26 


+ 66 + 1^3- 
1 + 26 2 + 36 


|666 + 


3?3 




-26 


• 66 • 


|666 + 


3^3 



A + 
A 2 . 



To^usc these identities to write_an expression for P' in closed form, intro- 
duce I: H VT*M defined by £(U) = i~ x l{U), and endomorphisms B* and B 
of VT*M and iJ given by B*{i~H) = r*V ' t - H s and B(U) = t^Vus. If A is 
taken to be cither (e/m 2 )B* or (e/m 2 )B then, since both maps have the same 
invariants, P'' L can be seen to be of the form 

P' ± = 

(l+Pl +P2^ 1 B*+p 3 (^B*) 2 )P 1 - - !H-( qi+ q 2 ^B* + q 3 (^B*) 2 )lP 

+4^(91 + q 2 ^B + g3 (J^B) 2 )r 1 P ± - (pi + P2 ^B* + P3 (^B*) 2 )P 
m z m z m z m z m z 

where pi, p 2 , p 3 , 51, 52, <Z3 have explicit rational expressions in terms of the 
invariants and have the following orders in e; p\ = 0(e 4 ), p 2 — 0(e 3 ), p 3 = 
1 + 0(e 2 ), gi = 0(e 3 ), q 2 = 1 + 0(e 2 ), g3 - 0(e). 

The Gauss equation states that if K s is the intrinsic sectional curvature of 
s(N), then K s is related to the ambient sectional curvature by 

K(s*U, s*W) - g(a(s*U, s*W), a(s*U, s*W)) + g(a(s*U, s*U),a(s*W, s*W)). 

To isolate the terms of order less than e that are present in this expression first 
decompose V s ,uS*W as V s ,uS*W — T(s*U, s*W) + S(s*U, s*W) where T is 
vertical and S is horizontal. Since V,s T = 0, it follows that B*i~ 1 s = 0, and 
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so the only terms of order less than e that contribute to K s arc of the following 
type 

g(a(s*U,s*W),a(s*V,s*Z)) = 

g(T(s*U, s*W), T(s*V, s*Z) + g{T(s*U, s,W), B*£S(s*V, s*Z)) + 
g(B*£S(s*U, s,W), T(s*U, s*Z) + 0(e) 

4 

777 

= -— g (U,W)g(V,Z) + 0(e). 

Hence, the Gauss equation gives that K s (s*U, s*TU) = —(m 2 /e)K(U, W) + 0(e). 
In fact to order e, a calculation shows that 

K s (s*U,S*W) = 

-—K(U, W) + ^{g(R(U, W)s, R(U, W)s) - g(V w R(U, W)s, V uS ) + 
g(VuR{U, W)s, V w s) + 3g{R{U, W)V uSl V w s) + 
ff(Viv(Vt/a), Vu(Vwa)) ~ g{V w (Vws), V^V^))) + o(e). 

Thus the terms that must be added to K(U, W) to obtain K s (s*U, s*W) occur at 
two order of e higher than K{U, W) and involve V.s, V.(V.s) and V.(V.(V.s)) 
as R(U, W)s = Vu(V w s) - V w (V v s). 

It is interesting to note that if GK S is truncated at order e 2 , and is used to 
construct approximate solutions for Friedmann geometries as in Section 2, then 
the solution curve Jo does not return to the infationary state but rather recedes 
to infinity in backward finite time for all values of a. 
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